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We study a theoretical model for the magnetothermal conductivity of a spin-| ladder with low 
exchange coupling (J <C @d) subject to a strong magnetic field B. Our theory for the thermal 
transport accounts for the contribution of spinons coupled to lattice phonon modes in the one- 
dimensional lattice. We employ a mapping of the ladder Hamiltonian onto an XXZ spin-chain in 
a weaker effective field B e fj = B — Bo, where Bo — - Bcl +- B <= 2 corresponds to half-filling of the 
spinon band. This provides a low-energy theory for the spinon excitations and their coupling to 
the phonons. The coupling of acoustic longitudinal phonons to spinons gives rise to hybridization 
of spinons and phonons, and provides an enhanced B-dependant scattering of phonons on spinons. 
Using a memory matrix approach, we show that the interplay between several scattering mechanisms, 
namely: umklapp, disorder and phonon-spinon collisions, dominates the relaxation of heat current. 
This yields magnetothermal effects that are qualitatively consistent with the thermal conductivity 
measurements in the spin-j ladder compound Br4(CsHi2N)2 (BPCB). 



I. INTRODUCTION AND PRINCIPAL 
RESULTS 

Quasi one dimensional (ID) magnetic systems are present 
in a variety of new compounds with magnetic elements, 
and provide interesting manifestations of strongly correlated 
physics in electronic systems^. These systems are realized in 
crystals with a chain-like structure of the magnetic atoms, 
where intrachain exchange interactions are much stronger 
than interchain interactions. Their low dimensionality leads 
to the enhancement of quantum fluctuations, and the forma- 
tion of exotic phases at low temperatures. 

In particular, spin- ^ chain systems (most commonly re- 
alized in Cu-based compounds^ are typically insulators in 
which the charge degree of freedom is frozen, and the dynam- 
ics is restricted to the spin sector. The elementary excitations 
are spin flips propagating along the chains direction. These 
can be described in terms of interacting Fermionic degrees 
of freedom, called spinons, which carry spin but no charge^. 
These systems therefore provide one of the simplest realiza- 
tions of Luttinger liquids (LL). This spinon LL is, in fact, 
the most abundant form of the so-called "spin-liquid" state, 
characterized by a magnetically disordered ground-state and 
power law spin-spin correlations^. 

The most elementary model for ID spin systems is the XXZ 
HamiltoniarP, describing a spin- 1 chain with nearest neighbor 
interactions, 

Hxxz — '^^Jxy(S? +1 S?+Sf +1 Si)+J z '^ i Si +1 Si—B'^2Si ■ 

i i i 

(1) 

Here J a > corresponds to antiferromagnetic exchange in- 
teraction, and B is an external magnetic field (note that here 
and throughout this paper we adopt units where g(j,B — = 
fi = 1). The isotropic case J xy = Jz yields the ID Heisen- 
berg model. On each site the spin operator is represented by 
Si = where a a are the pauli matrices. The spin chain can 
be mapped into interacting spinless Fermions on a lattice^ 121 , 
where the magnetic field B serves as a chemical potential. At 
zero field the Fermions are at half filling, and upon raising 
the magnetic field they gradually polarize until saturation at 
B — B c which corresponds to a depletion of the spinon band. 

More complicated variants of the XXZ and Heisenberg 
model can describe quasi ID systems with additional in- 



teractions such as zig-zag chains, spin-Peierls chains, and 
ladder!®^. These systems support a richer phase diagram 
including, e.g, gapped dimer crystal phases. Upon tuning 
the magnetic field the system may undergo a phase transition 
from a gapped phase into a spin-liquicP. In particular, in a 
ladder subject to a strong field B c \ < B < B C 2, a LL phase 
of gapless spinons is recovered^. 

One of the prominent manifestations of a spin-liquid state 
is the contribution of gapless spinons to transport. Since there 
is no straightforward way to measure the spin current through 
an antiferromagnetic chain, investigation of the spinons prop- 
erties can be done by measuring the thermal conductivity k. 
Experimental evidence for a substantial enhancement of ther- 
mal conductivity along the chains direction (km), has indeed 
been found in CuO based chain compound s 10 " 12 !. However, 
interpretation of the data is complicated by the dominant 
contribution of crys tal ph onons, and in particular their cou- 
pling to the spinoniP^3] In principle, an obvious means of 
disentangling the spin degrees of freedom is the application 
of an external magnetic field B, which allows the tuning of 
system parameters in the spin sector only. The resulting mag- 
netothermal effects - namely, variations of k as a function of 
B - can serve as a valuable probe of the spin system. At 
low temperatures both spinons and phonons contribute to the 
heat transport. The total heat conductivity can be split into 
a pure phononic contribution, K p h(T), and a magnetic part, 
Kmag(B,T). Then, we can extract the magnetothermal con- 
ductivity: 

Ak(B,T) = k(B,T) - k(0,T) = K mag {B,T) - K ma9 (0,T) . 

(2) 

Magnetothermal effects as mentioned above are practically 
inaccessible in the typical CuO compounds, where the large 
exchange coupling (J of order 2000 K) dictates an enormous 
scale of the desired external field. In contrast, a field-tuned 
manipulation is easily accessible in organic based magnetic 
compounds, where J is typically of order 1QK. An experiment 
in the organic spin- chain material Cu^tLNaXNOs)^ mea- 
sured the magnetothermal conductivity. It indicated a non- 
monotonic B-dependence of km, and in particular - a pro- 
nounced dip feature with a minimum at a field scale B m in ~ 
T. A subsequent theoretical study^has shown that such fea- 
ture arises due to the interplay between disorder and umklapp 
scattering of the spinons: the latter process is sensitive to the 
field-induced tuning of the spinon Fermi-level away from the 
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middle of the band. It thus reflects the Fermionic character 
of the spinons. 

As opposed to the spin-chain compounds mentioned above, 
in spin- ladder compounds, magnetothermal effects are ex- 
pected to dominate at high B where the spin-gap closes up. 
A recent experiment!^ measured the magnetothermal conduc- 
tivity in the spin-ladder compound Br4(CsHi2N)2 (BPCB). 
An experimental study^ of thermodynamic properties of this 
compound confirmed that it is described very well by the spin- 
ladder model with J« = 3.6K (the exchange along the legs of 
the ladder) and J± ~ 13K (the exchange along the rungs), 
and its appropriate LL representation in the gapless regime 
B c i < B < B c 2- 

Indeed, the experimental data of Ref. [17] indicate that 
upon raising the magnetic field, the magnetothermal conduc- 
tivity Ak(B) vanishes for fields smaller than B c \. However, 
when the magnetic field is raised further and the spin-gap is 
closed, there is a large decrease in the magnetothermal con- 
ductivity. On top of this decrease there is a double dip fea- 
ture with a local maximum at Bo = Bcl + Bc2 , corresponding 
to half-filling for the Fermionic excitations. We assert that 
this data can be qualitatively explained as follows: first, the 
spinons in this system are slower than the phonons, there- 
fore they act as impurities for the phononpSI. This effect 
induces a decrease in the conductivity upon entering the spin- 
liquid regime (B c i < B < B^). Second, around half-filling 
(B = So) there is a positive spinonic contribution to trans- 
port observed as a maximum at B — Bo . The doub le dip fea- 
ture resembles results obtained for spin- chain i 15 ^ 16 ^ where the 
minimum in the magnetothermal conductivity corresponds to 
moving the chemical potential away from half filling, to a scale 
of order T. 

Motivated by these observations, in the present paper we 
study a minimal model for the magnetothermal transport of a 
coupled spinon-phonon system in a single ladder. Our theory 
accounts for a crucial distinction between ladders and chains: 
the strong magnetic field required to enter the gapless spinons 
phase provokes an enhanced coupling between spinons and 
phonons. This leads to hybridization between the spinons and 
phonons excitations. In addition, scattering of the phonons by 
the slower spinons is magnified, generating a relatively strong 
negative contribution to Ak(B). Qualitatively, our calculated 
Ak(B, T) resembles the experimental data of Ref. [17j . 

The paper is organized as follows: in Sec. [II] we derive 
the low-energy model for the spin system in the presence of 
coupling to ID phonons. In Sec. |III| we study the effect of 
scattering processes on the thermal conductivity in the frame- 
work of the memory matrix approach for the calculation of 
the conductivity tensor, and obtain the leading magnetic field 
and temperature dependencies of the thermal conductivity k. 
In Sec. |IV| we summarize and discuss the results. Finally, in 
appendices A through C we present details of the calculation 
of the various memory matrix elements. 



where t = J^ f /2 and V = J e J s . For B eff = the 
Fermionic band is half-filled and the Fermi momentum is 
= jjL. Finite B e ff corresponds to a chemical poten- 



II. LOW-ENERGY MODEL FOR THE 
COUPLED SPIN-PHONON SYSTEM 

We wish to compute the thermal conductivity of a system 
which consists of antiferromagnetic spin-| ladders interacting 
with the lattice phonons. To this end, we focus on a simplified 
model for such a system, which considers a single ladder - i.e., 
both spinons and phonons are one-dimensional. The param- 
eters of the model are adjusted to mimic those of 
in particular assuming the limit J± > Ju (strong rung cou- 
pling). In addition, we assume Jy <C Qd- In this section 
we describe the low energy model of the system, and derive 
the eigenmodes which constitute the elementary excitations 
of the coupled spin-phonon system. 



A. Bosonization of the spin ladder Hamiltonian 

We begin by describing the spin system. The Hamiltonian 
of a spin- 1 two leg ladder in a magnetic field B along the 
2-direction is 

JV 2 N 

H a = y~^[J||Sj,i/ ■ Sj+i,!/ — BSl tV ] + J±Si,i ■ Si, 2 , (3) 

■i— 1 u=l i — 1 

where v = 1, 2 denotes the leg index. For J± > Ju, it can 
be approximately mapped into an effective spin-| chain in a 
weaker magnetic field™ B e ff = B — Bo: 

N N 

(4) 

where the effective parameters are given by 

Jxy Jz ^11 ' 

B eff = B -B , B = Jx + J||/2. (5) 

The isospin operators erf describe the effective spin-| dynam- 
ics characterizing the low energy sector, which at high B is 
restricted to the singlet and lower triplet state on each rung. 
Hence, in distinction from the real-spin XXZ model [Eq. , 
(of) = corresponds to a time-reversal symmetry broken 
state. To derive the low-energy model for the dynamics of 
this system we first use the Jordan- Wigner transformation, 

i-l 

<7+ -> ct exp(i?r c ] c i)i °i -> c ! c * _ V 2 (6) 

j = -co 

which maps the spin problem onto a model of interacting 
spinless Fermions on a lattice: 



(7) 



tial for the Fermions, which shifts the Fermi momentum into 
fc F = fep (1 + M e ff), with M e ff an effective magnetization. 



H s = -tJ2(4^+i + h.c) +Vj2(t$Ci - l/2)(ct +1 c i+ i - 1/2) , 

i i 

I 
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Near the middle of the band (Be// = 0), the Fermion oper- 
ators can be expressed in terms of Bosonic ones related to the 
Fermion density fluctuations using the standard dictionary of 
abelian Bosonization (see, e.g., appendix D in Ref. [I]). For 
the spin operators (in the continuum limit: x — ia) this yields 



iB( X ) 



a+(x) = 



a z (x) = d x <j){x) + 

1 + 



-[(-)- +cob(20(x))] 



l 2 



cos(2<j!>(:r)) , 



(8) 



where o ± (x) = 775 c*, ° z {. x ) — , an d a is the lattice 

constant. Substituting Eq. (JsJ) into Eq. Q we can describe 
the low energy properties of the spin system in terms of the 
Boson Hamiltonian: 



H a — H g + H u , 
1 

2^ 



HI 
H u = g u 



dx[g(d x (f>(x)) 2 + v F (nU(x)) 2 
dx cos[40(a;)] , 



(9) 



where 



v F = aJ\\ , g = v F ( 1 + - ) , g u = 



in 



4n 2 a 



(10) 



and II(x) = -9^0(3! ) is the canonical conjugate of 4>(x), obey- 
ing [n(a;), 4>{x')\ = iS(x — x'). H® is the standard LL Hamil- 
tonian 



Hll = 



2tt 



dx 



K 



{d^f + K(nliy 



(11) 



where 



v F [l+ - 
n 



has the dimensions of velocity and 



K 



1/2 



-1/2 



is the dimensionless Luttinger parameter. Since K > 1/2, 
the umklapp term H u is irrelevant (i.e. flows to zero under 
renormalization group (RG) for T — > 0) and hence can be 
neglected in the description of the low-energy thermodynamic 
properties. However, as we shall see in the next section, it 
plays an essential role in the transport. 

For B 7^ Bo, the finite Be// introduces an additional term 
to due to the last term in Eq. ([4]), which induces a finite 
effective magnetization. The most relevant correction is of 
the form 

i J dxB eff d x <j> , (12) 

which can be absorbed in the Gaussian part by a shift of the 
field cf>, reflecting the shift of chemical potential for spinons. 



As implied by the exact Bethe ansatz solution, the LL form 
of Hg is in any case maintained for arbitrarily large B e / / , 
but with renormalized parameters'^. In particular K(B e ff) 
approaches 1 close to the edges of the band (B — > B c i or 
B -> B c2 ). 

An additional correction to H s arises from weak disorder in 
the lattice, which can be accounted for by adding a random 
term 5B(x) to B e ff in Eq. Q. Such term may arise from 
defects leading to random corrections to Jj_, J« via Bo [Eq. 
d5j]. This introduces a scattering term proportional to 



dxSB(x) cos[2<?!>(:r)] 



(13) 



When we discuss the relaxation of the heat current, both the 
umklapp and disorder terms will become important, and will 
be considered as perturbations of H® . 



B. Coupling to Lattice Phonon Modes 

Up to now we described only the spin system. Next we will 
include the phonons in the model. In a single two-leg ladder 
of atoms, three modes of ID phonons should be accounted 
for: two acoustic modes, longitudinal and transverse, and an 
optical mode associated with fluctuations in the rung length. 
The dominant coupling of phonons to spinons arises from the 
dynamical corrections to the exchange interaction, J, due to 
lattice vibrations. The spinons therefore couple to leading 
order only to the longitudinal acoustic mode (via fluctuations 
in J||) and to the optical mode (via fluctuations in J±). 

We first consider the effect of coupling of optical transverse 
phonon modes to the spinons. We show that such coupling 
merely leads t o no rmalization of the Luttinger parameters, u 
and K of Eq. |TT|. 



Let us define the transverse phonon field in the following 
way: 

Ut(x) = U 1 (x)-U 2 (x), (14) 

where Ui : 2 are transverse displacements (along the rung di- 
rection) of atoms in different legs of a ladder normalized by 
the rung size 6. Now, we substitute this definition of U t in 
the phonon-dependant exchange to obtain: 



J±(U t ) « Jl >+5J ± U t , SJ± 
r = b(l + U t ). 



. dJ± 
dr 



(15) 



Inserting into Eqs. Q and ( 12 1, we find that this adds to the 
Hamiltonian a term of the form 



H sp = —5J± 



Ut(x)d x 4>(x). 



(16) 



It is useful to change into momentum representation with 



Ut(x) = J J ^> fc e fc + b\e-^), fa) = ^> fc e lte + ^e^) (17) 



k 



r 



(where L is the length of the legs). Then, the quadratic part of the coupled spinon-phonon Hamiltonian is given by 

H — H° + Bp + HI p , (18) 
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where H° describes the spinons in terms of a Luttinger Hamil- 
tonian [Eq. (Ill], 



(19) 



describes the optical transverse phonons, and 

Hi p = -5J±aY J iK<t>kb\-<t>-kbk) (20) 

k 

is the spinon-phonon interaction. Using a coherent path in- 
tegral representation, it is a straightforward exercise to inte- 
grate over the phonon degrees of freedom, yield an effective 
action for the spinons, S e ff, defined as 



= / DbDbe 



-S[<t>,b,b] 



In the limit k, uj n — > this results in a Luttinger model with 
a modified coefficient of (d x <f>) 2 : 



u 
K 



u SJ, a 
K cj 



Hence, the renormalized parameters become: 



K — K 



1 



&j\aK 



1 -1/2 



1 + 



SJlaK 

LJQU 



-,1/2 



(21) 



(22) 



Next we focus on the longitudinal phonons, which coupling 
to the spin sector has the most dramatic consequences. As- 
suming small displacements of atoms from their equilibrium 
positions, we can approximate the exchange interaction by: 



J|,(r)wJ,j 0) +g ]l d x U 1 (x), si. 



a 2 ^\ 
a Qr \ r=a , 



(23) 



where r — a[l + Ui(x + a) — Ui(x)] is the distance between 
neighboring atoms on the same leg, and the dimensionless 
field Ui(x) describes the relative longitudinal displacements of 
atoms. When inserted into Eq. Q, these corrections give rise 
to coupling between the spinons and phonons. The Hamilto- 
nian describing longitudinal phonons traveling parallel to the 
chains is 



v 

2tt 



dx[(nPi(x)Y + (d x Ui(x)) 2 



(24) 



where v ~ aQjj (with 9d the Debye temperature) is the 
sound velocity, and Pi is the momentum conjugate to Ui. 

After inserting the phonon-dependant corr ection to the ex- 
change interaction into Eqs . |5| , (JoJ) and ( 12 I, and adding the 
phonon Hamiltonian [Eq. | |24| ] , the low-energy Hamiltonian 
of the coupled spin-phonon system can be written as: 



H = 



1 

2^ 



dx {g(d x 4>(x)) 2 + v f (ttII(x))' 2 + hid x Ui(x)d x (j>{x) 



B (o) 

+v[(dM(x)) 2 + (itPi(x)) 2 } + y ±^Hd x Ui{x) 



(25) 



with 
hi = 



-9\\ 



1 + 



2(1 + 2/ir)B^ f a 



d x <j){x) = d x <j>(x) + 



B 



(0) 



B 



(0) 
eft 



B-J 



(o) 



j{ 0) /2. (26) 



In Eq. (25 I we neglected small terms (of order d x Ui(d x (j>) 2 
and higher) . These terms are irrelevant, and moreover corre- 
spond to forward scattering that cannot contribute to trans- 
port properties of the spinons to leading order. Note that, 
in contrast with spin- chainsP 1 , at half-filling (B ( ° f ] f = 0) the 
coupling to the phonons via the coupling constant hi is lin- 
ear in the spinon field d x <f> and has to be included in the 
low-energy Hamiltonian. This reflects the breaking of time- 
reversal symmetry in the system, where spinons correspond 
to fluctuations around a partially polarized magnetic state. 
Below we show how these terms lead to new eigenmodes of 
mixed spinon-phonon degrees of freedom. 



C. Derivation of Hybrid Eigenmodes 



The Hamiltonian Ho in Eq. |25| describes the low energy 
properties of the coupled spinon-phonon system. In order to 
find the eigenmodes which constitute the elementary degrees 
of freedom of the system, we proceed in diagonalizing it by a 
canonical transformatiorP^l 



4>{x) = C4>i(x) - X 'Sfcix), 

Ui(x) = —SfaW + Cfcix), (27) 



and similarly for the canonically conjugate momentum: 



where: 



n(a;) = CTIi(a;) - ^Sn 2 (x), 

Pi(x) = A 2 OT 1 (a;) + C , n2(x), (28) 



C ■: 



1 

71 



A 2 = 



1 1 1/2 
v / 3 5 TlJ ' 

hi\Jv F /v 
gvF/v — v 



A ■ 



1 

= 71 


i + 






v V 


V 



1 



-,1/2 



< 1. 



(29) 



Eqs. (|27j) and ([28]) are designed to preserve the canonical 
commutation relations [<j> v {x) , H v >(x')] = i8„„i5(x — x'). The 
last approximation in Eq. ( 29 1 assumes vp v, which follows 



from J|| <C @ d ■ The parameter A defines the strength of 
the coupling between the spinons and phonons. Note that it 
would be much stronger in a compound where J ~ Qd, in 
which case the phonon and spinon velocities match, vf ~ v. 



After this transformation, Ho [Eq. ( 25 1] takes the form 



J Lf=l,2 L 



+ - 



with 



hiB 



d x (j>i(x) + Cd x 4> 2 {x) 



(30) 
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Vl 1 

Ki ~ 2V A 2 + 1 

V2 __ 1 
K2 ~ 2\M 2 + 1 



l(VT+A~2 - 1) + — + 1) + .P^-hiA 



gvF 

V 



(Vl+A 2 + 1) + «(Vl + A 2 - 1) - hnl—A 



, U 2 #2 = «, 



(31) 



which can be approximated for vf ^ v by 

vi ^ v 2 L SA 2 \ v 2 ^ gv F L 3A 2 v 2 
K\ vf \ 4 / ' iT 2 w V 4g«F 

In this form, _ffo is separable into two independent species of 
LLs. Using A <C 1, the LL parameters are approximated by 



(32) 



v 1 



3A 2 
4 



1/2 



1 + 



1 + 



4gvF 



vf 
v 



1 + 



3A< 



1/2 



A" 2 



'I'F 



\ 



1 + 



igv F 



1 + 



(33) 



Finally, to get rid of the linear terms d x <j>i( x ) an( A d x (f>2(x) in 
Eq. |30l , we define 



7/\ , / \ , B e ffKlhlS 

M*) = M*)+ 2vigX2 * 

B e ffK2hiC 



4>2{x) = <fa(x) 



2v 2 g 



and 



ni(as) = ni(a:), U 2 (x) = U 2 (at) 



(34) 



(35) 



which preserve the canonical commutation relations. The low 
energy Hamiltonian is now cast in the quadratic form of a LL: 



2 r 

H = Y. V » / 

i/=i 



dx 
2^ 



K, 



-{d x <t> v {x)f + K„(*n„(x)y 



(36) 



The Hamiltonian (361 is integrable (i.e. it has an infinite 
number of conservation laws), therefore the currents we are 
interested in (e.g heat current) are protected and cannot de- 
grade. In order to get a finite conductivity we must add per- 
turbations around Ho, e.g. the previously neglected umklapp 
term, which in terms of the shifted spinon field is given by 



H u = g u I dx cos[4<j>(x) — Akx], Ak 



4B t 



(37) 



H u describes processes where two spinons move from the right 
Fermi surface to the left (or vice versa), gathering momentum 
Afc = 4&f — G in which G = ^ is the reciprocal lattice mo- 
mentum. Another important correction to Hq is the backscat- 
tering term 



H d 



dx£(x) cos[2^>(ie)], <^(x) 



SB {2 



(38) 



which describes scattering of spinons due to weak disorder 
caused by defects in the lattice. We assume uncorrelated ran- 
dom disorder where the sample average gives 



C(aO = 0, 

((x)((x') = D5(x - x') . 



(39) 



Using Eq. (27 1, H u and Hd can be expressed in terms of the 



hybrid spinon-phonon eigenmodes </>i(x), <p2{x) 

H u = g u J dx cos(2a</>i (a;) — 2/302 (x) — Akx), 
H d = [ dxC(x)cos(a<j>i(x)- P</> 2 (x)), (40) 



with 



2C, 13 = 2\ 2 S . 



(41) 



Finally, we note that higher orders in the expansions Eq. 
( |15[ ) and ( |23[ ) yield an additional scattering term between 
phonons and spinons, which turns out to have a significant 
effect on the transport carried by phonons. subsection |III B| 
is devoted to a detailed study of the implication of this term 
on the th erm al conductivity. Together with Hd, H u [Eqs. 
( 37 1 and 1 38 \] , this scattering process governs the degrading 
of currents leading to a finite conductivity. 



III. THERMAL TRANSPORT OF THE 
SPIN-LADDER SYSTEM 

The magnetothermal effects observed, e.g., in Ref. [IT] 
are a consequence of the interplay between different scatter- 
ing mechanisms which result in the change of the thermal 
conductivity as a function of magnetic field [n(B)]. Two pri- 
mary effects are expected to dominate the B-dependance in 
the coupled phonon-spinon system: one arises from the pos- 
itive contribution of spinons as heat carriers, and the other 
from the negative contribution of spinons acting as scatterers 
of the phonons. The former contribution is governed by the 
interplay of two scattering mechanisms, umklapp and disor- 
der, and consequently depends on the deviation of the spinon 
chemical potential from a commensurate valued At the same 
time, the scattering due to phonon-spinon interaction is also 
dependant on the filling of the spinon band, which dictates 
the available phase-space for scattering. We note that due to 
the hybridization of phonons and spinons, these various effects 
are not entirely separable. In the present section we derive 
the magnetothermal conductivity of the spin-ladder model by 
using a memory matrix formalism, which allows an account of 
all the above mentioned scattering processes on equal footing. 



Approximate Conservation Laws and Memory 
Matrix Formalism 



To calculate the heat conductivity of the spin-ladder sys- 
tem, we use the memory matrix formalisnP^ which has been 



G 



successfully imple mented in previous studies of thermal trans- 
port in spin-chain s 13 ^ 1 . The memory-matrix approach is 
suited for systems where due to approximate conservation 
laws, the conductivity almost diverges^. The main step 
within this approach is the calculation of a matrix of re- 
laxation rates for a given set of slow modes. The method 
allows to calculate transport coefficients within a hydrody- 
namic approximation, and provides a reliable lower bound to 
the conductivity^. Moreover, it gives precise results as long 
as all the relevant slow modes are included in the calculation. 

The heat transport properties of the spin-ladder system at 
low temperature are governed by the approximate conserva- 
tion of a certain current, J c which has a finite overlap with 
the heat current. In particular, an exponentially slow decay 
of J c will lead to an exponentially large heat conductivity: 
the component of the heat current overlapping with J c is pro- 
tected, and will decay exponentially slowly. 

In the present case, the important step is to realize that in 
the absence of disorder, the linear combination 



Jc — J d) 



Ak 



Js 



(42) 



is conserved, meaning [J c ,Ho + H u ] — 0, where Ho is the 
LL Hamiltonian (361, and H u the umklapp term [Eq. ( 3T l] . 
Here, 



J,i, 



dxll(x)d x (j)(x) 



J S = N R -N L = j dxU{x) 



(43) 



are the (normalized) heat current associated with the spinons 
and the spin current, respectively, where Nr and Nl are the 
total number of right or left moving spinons. The overlap 
between the heat current and the conserved current is man- 
ifested by the appearance of J$ in Eq. (42 1. The reason 



that J c is conserved by umklapp scattering is as follows: the 
umklapp term describes a process where a momentum Ak is 
generated and therefore induces a change in J$ proportional 
to Ak. In the same process, the normalized spin current is 
changed by —4 as two right-moving spinons are scattered into 
left-moving states. Since J c is conserved by the umklapp term 
H u , the heat current can not be degraded by H u alone and 
additional scattering processes need to be accounted for. 

The low-energy Hamiltonian Ho conserves an infinite num- 
ber of modes in addition to J c . However, when perturbations 
are added, these modes decay faster than the conserved cur- 
rent Jc, since these modes do not commute with all the terms 
added. 

We now show how to calculate perturbatively the thermal 
conductivity when the relaxation of the heat current is dom- 
inated by the slow modes. In our case the memory matrix is 
formulated in a space spanned by the slow modes Ji, J2, J3 
and J s [Eq. J43b], where 



Ji = j dxlli{x)d x <f)i(x), J 2 = j dxU 2 (x)d x (j}2(x), 
Js = / dxn 3 {x)d x Mx) (44) 



which are all conserved by Ho- The fields <f)3^3 represent 
the transverse acoustic phonons which do not hybridize with 
spinons but are still scattered by spinons and therefore con- 
tribute to relaxation of the heat current. The heat current 



along the chains direction is Jq = v\ Ji + v\j 2 + ff J3, where 
vi and V2 are given in Eq. ( 33 1. 

To set up the memory matrix formalism^} we first intro- 
duce a scalar product on the operators in the space spanned 
by the slow modes 

(A{t)\B) = T f 1/T d\(A 1 (t)B(i\)), (45) 
Jo 

where (...) denotes an expectation value at equilibrium, in- 
cluding average over disorder configurations. The dynamic 
correlation function of the operators A and B is 



iT 



dte^\A(t)\B) 



dte lLJt < [A(t),B] > - 



(A\B) 



and the matrix of conductivities is given by 



(46) 



(47) 



where p, q are either of the slow modes. The heat conductivity 
is given by 

1 



(48) 



where Jq denotes the heat current. One can also write the 
matrix of static susceptibilities as 



Xpq 



TL 



(Jp\Jq 



(49) 



It can be showrP that the matrix of conductivities, a, can 
be expressed in terms of a memory matrix, M 



a = x(T)[M(u,T) -iuxiTT'xiT) 



(50) 



The elements of the matrix M{uj) in the d.c. limit (u — > 0) 
are 



M va = lim 



C pq (oj) - C pq (u = 0) 



-id u C pq \ w=0 (51) 



where p, q can be each of the slow modes of the theory, and 
C pq (uj) is the Fourier transform of the retarded correlation 
function, 



C pq {io) 



dte-\[F p {t),F q {Q)\), 



of the force operators: 

F p ~ j p = i[H,J p ]=i[H pert ,J p ]. 



(52) 



(53) 



Here H per stands for perturbations to the low-energy Hamil- 
tonian of the system Ho, that can relax the current J p [such 
as H u , H d , Eqs. ( 37 1 , (38 (]. In the last equality we used 
[J p ,Ho] = for p = s, 1,2,3, which justifies a perturbative 
expansion of M: since J p are already linear in perturbations 
around Ho, the expectation values in Eq. (49 1 and (521 are 
computed with respect to Ho- 



conductivity is given by 



From Eqs. (481 and (50 1 it follows that the d.c. thermal 

(54) 
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The static susceptibility matrix is given by 









X = 



V 



ttT a 









o 



(55) 



(where the matrix indices are s, 1, 2, 3). In our case the three 
leading perturbations contributing to the memory matrix, M, 
are umklapp and disorder in the spin sector, and phonon 
scattering processes which include phonon-spinon interaction. 
Thus, the memory matrix is separated into three parts, 



(56) 



M = M u + M d + Ab- 



using the conservation law of the slow mode J c [Eq. |42|] we 
find relations between the different umklapp matrix elements 
(see appendix A): 



Aft = 


M" 

4 ss 


-AfS, 




Af? a = 


Ak ,« 


- M 2 " 2 , 


(57) 


Mu = 




- MS- 





When < «, and thus 7\ia 2 « if 2 /3 2 [see Eq. j4l|l, the 
matrix M° greatly simplifies. Using the relations (571 we 
find that the leading contribution to M u depends only on one 
element M?, = Mf. and we have 



M u ^ 



Ml 




4 



A/, 



\ 





^M s " (^) 2 M S " 



\ 



(58) 



0/ 



The disorder contribution, M , is a diagonal matrix with 
elements denoted by M d p = M d (p — s, 1,2, 3). Finally, 
the dominant contribution from phonon-phonon and phonon- 
spinon scattering, Al p , appears in the diagonal elements Mi, 
Mz (a detailed calculation is provided in the next subsection) . 



Substituting these relations into Eq. 
expression for the thermal conductivity 



54b. we obtain an 



k(B,T) 



1 

T 



"42 

Mi 



+ 



164xl 2 (M a d + MD 



Afc 2 M s d M« + 16M$(M$ + M s " 



+ 



4 2" 

"3X33 

M, 



(59) 



The B-dependance of this expression is encoded in the vari- 
ous matrix elements M", Mf , Mi and M 3 . We note that the 
scattering processes in the spin sector (described by M™, Mf) 
dominate near half-filling of the spinon band, where their 



spectrum can be linearized and Bosonization is justified. Us- 
ing Eq. ( 40 I for the relevant terms in the Hamiltonian, we 



derive M" (for a detailed calculation, see appendix A), 



/ rp \ 2K 2 (S 2 +2K ia 2 -3 

M: =9s{y) B[K 2 p 2 /2 K 2 P 2 ]B[K 2 P 2 /2 + iS,l- K 2 p 2 ]M[t>(l - K 2 p 2 /2 - iS) - {<S>{K 2 p 2 /2 + iS)] 



_ gi ■ i ts o2n <■ _ v 2 Ak _ v 2 



(60) 



r 



where B(x,y) — ^^2+^ * s ^ ne Beta function, ^(x) — j4^j 
is the digamma function and the parameters K%, K 2 , a and 
/3 are defined in Eqs. (1331) and 141 h - The dimensionless pa- 



rameter 8 deter mine s the dominant field dependance of M" 
via Afe [see Eq. ( 3T l] - For the disorder part of M we find 



K 2 /3 2 / 2+ K ia 2 / 2-2 





K 2 p 2 /2+K ia 2 / 2 






:«) 





K 2 2 /2+K ia 2 /2 



o\ 




0/ 



(61) 
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with 



D a S 



4nDa 2 



2t;i 



Do S 



2d 2 



(62) 



where D is defined in Eq. (391. The B-dependance of M is 
implicit in the parameters K\, K2, a, /3. 

Incorporating Eqs. ( |60[ ), ( |61[ ) in the second term of Eq. 
( 59 1 , we find the primary magnetothermal effects originating 



from relaxation processes in the spin sector. To complete the 
derivation of k(B, T), one must account for the B-dependance 
of Mi , M3 originating from phonon-spinon interaction. This 
part of the derivation requires special attention, and is dis- 
cussed in subsection B below. 



B. Phonon-Spinon Scattering 

As already mentioned above, we focus on ladders with small 
exchange coupling obeying Jy «6d, where the phonons typ- 
ical velocity is much larger than the spinons velocity. There- 
fore, the spinons act as impurities which scatter the phonons. 
These scattering processes lead to relaxation of the phononic 
heat current and therefore to a prominent dip in the thermal 
conductivity upon entering the partially filled spinon band, 
for B c i < B < Bc2- These add a B-dependant contribu- 
tion to the scattering rate of both longitudinal and transverse 
branches of acoustic phonons (represented by the Bosonic 
fields Ui,(j>3). Below we study the contribution of phonon- 
spinon scattering processes to the corresponding memory ma- 
trix elements. 

In order to account for scattering processes of phonons on 
spinons, we expand the phonon-dependant exchange ( 23 1 to 
second order in d x Ui. 



The most relevant phonon scattering 
term arising from this order of the expansion is of the form 



Hps = gp-s J -^(d x Ui) 2 d x 



(63) 



A similar term resulting from the transverse displacements 
yields a coupling of spinons density to (dxfis) 2 - We wish to 



consider magnetic-field dependance in a wide range, i.e. the 
entire spinon-band. For this purpose we model the spinons 
as free Fermions at a chemical potential dictated by B, a rea- 
sonable approximation, e.g., for the ladders in BPCB where 
the Luttinger parameter K is not far from 1 for arbitrary B 
[9]. Using d x (j)(x) = — Trp(x) pQ, and turning to Fourier space, 
H P s acquires the form: 

Hps = ^2 V pp 'c\ l c k b [ v ,b p 8{k - k' +p-p -G n ), 

kk' pp' G n 



V m 



sPP 



yJOJpLUpl 



(64) 



where c k is a Fermionic (spinon) annihilation operator, and b p 
is a Bosonic (phonon) annihilation operator. H P s describes 
elastic scattering where a phonon and a spinon with momenta 
p and k respectively, scatter into p' and k' respectively. 

We need to calculate the effect of this term on the phononic 
heat current 



Jph — ^ ' v p 



ujpbpbp, 



dujp 
dp 



usign(p), (65) 



where we have assumed a linear dispersion lu p « v\p\. The 
memory matrix element [Eq. ( 51 1] is therefore calculated with 



the correlation function ( 52 \ of the force operator 



F„ 



\_Jphi Hp 



(66) 



Using Wick's theorem we obtain an expression for this corre- 
lation function (see appendix C for details) 

C P s = ^2 W pq 5(Au))np +q fk- q (l + n v )(l — f k ), 



kpq 

Wpg = -2v 2 g 2 p_ 8 \p(p + q)\q 
Auj = e k — £k- q +u p — W. 



(67) 



+ l)" 1 , n p = (e^ /T - I)" 1 , are Fermi 



where f k = (e^-^/ T 
and Bose distributions respectively. The memory matrix ele- 
ment is the derivative of C p _ s with respect to uj. Integrating 
by parts and using the energies delta function we obtain: 



Mp 



I 



dkdp 



i')q 



W pq n p+q f k - q (l + n p )(l - fk) 



3q 



d q \1=1a 



(68) 



r 



where go is the momentum transfer that obeys the energy 
and momentum conservation dictated by the delta functions. 
Since the phonon dispersion ui p = v\p\ is much steeper than 
the spinon dispersion e k = — J|j ' cos(fca) , energy and momen- 
tum conservation can only be satisfied by phonon backscat- 
tering where p — > p + q — —p + Sp, \5p\ « \p\. This is 
because a small change in phonon momentum will lead to a 
large energy change, while the spinon energy transfer is small 
for small momentum transfer (see Fig. [l|. The integrals in 
Eq. ( 68 1 were solved numerically after approximating 



5p = — — [cos(fc + 2p) — cos(fc)] , 
v 



to get the temperature and field dependencies of M p _ s (see 
Figs. The memory matrix is closely related to the 

relaxation time of the scattering process M p _ s ~ i~ p ^ s . Indeed 
as indicated by Fig. [2] phonon scattering occurs practically 
only for B C 2 < B < B c \, in the spin-liquid phase where the 
spinons are gapless, and it is maximal for B — Bo (half filling 
of the spinon band). The temperature dependance of M p _ s 
(Fig. [3| for B e ff = gives a good fit to a power law M p - S ~ 
T a , with a ^ 4.5. 



We next recall that under the approximation vf «C v, one 
obtains a < f) [see Eqs. f29}, pT} 

while ( 



hybrid mode tj>i is phonon- 



ikc, 



which implies that one 
62 is spinon-like. There- 



fore, the scattering of longitudinal phonons is included to a 
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Momentum 



FIG. 1: (color online) Schematic Plot of the phonons (red) 
and the spinons (green) energy dispersions, the allowed 
phonon backscattering (red curved arrow) and the corre- 
sponding spinon scattering (green curved arrow). The black 
arrows represent the momentum transfer go = p' — p = k — k'. 



52 0.0004 - 

•a 

j3 0.0003 
1 0.0002 - 




B[T] 

FIG. 2: (color online) Isotherms of M p - S as a func- 
tion of magnetic field for various temperatures T = 
0.25K, 0.27K, 0.30K, 0.33K, 0.36K, 0.38K, 0A2K. The pa- 
rameters used for this plot: Jn = 3.67^, J± = 14.4K, g p - s = 
0.04, z = 18A". 



good approximation only in the Mi element of the memory 
matrix. This is added to the disorder term already retrieved 
earlier, and a B-independent contribution which assumes a 
power-law dependance on T. We thus obtain an expression 
of the form 



Mi = Di 



+ D P T'< + M v 



(69) 



A similar expression, excluding the first term, holds for M3 
which describes the scattering of transverse acoustic phonons. 



Substituting in ( 59 1, we obtain the final expression for k(B, T) 
and consequent y for Ak(B,T) [Eq. Q]. The resulting B 
and T dependance of Ak(B)/k(0) are plotted as a function 



'2 0.020 
3 



T[K] 

FIG. 3: (color online) M p _ s at half-filling (B e ff = 0) as a 
function of temperature. The parameters used for this plot: 
Jn = 3.6K, J± = UAK, g p - s = 0.04, v - = \%K. 



of magnetic field for different temperatures (Fig. |4|. We note 
that this result, although based on a highly simplified min- 
imal model which captures the main physics of the system, 
qualitatively reproduces the prominent features of the exper- 
imental data of Ref. 1171 



0.00 



o 



pq-0.15- 




B[T] 

FIG. 4: (color online) Isotherms of the normal- 
ized magnetothermal conductivity as a function 
of magnetic field for various temperatures T = 
0.25if, 0.27K, 0.30K, 0.33K, 0.36K, 0.38K, QA2K. The pa- 
rameters used for this plot: Jm = 3.6K, J± = 14AK, g p - s = 
0.04, D = 1.6 X 10" 8 , I = D p = 7.4 x 10" 8 , D z = 

4.6 x 10" 6 , 7 = 4. 



IV. SUMMARY AND DISCUSSION 

In this work we studied the thermal conductivity of weakly 
disordered spin ladders subject to a magnetic field and cou- 
pled to phonons. We found that due to coupling between 
the phonons and the spins, the elementary degrees of free- 
dom are hybrid spinon-phonon modes and strong scattering 
of phonons on spinons is induced. Our study of the phonon- 



10 



spinon scattering found that due to energy and momentum 
conservation only certain backscattering processes are allowed. 
The phonon-spinon scattering along with umklapp and disor- 
der scattering lead to a prominent dip in the thermal con- 
ductivity. We examined the mechanisms responsible for the 
relaxation of the heat current, and showed that an interplay 
between umklapp, weak disorder and phonon-spinon scatter- 
ings underlies the transport properties at low temperatures. 
For this system it leads to minima in the thermal conductiv- 
ity isotherms when the effective field is of the order of the 
temperature |-B e //| ~ T, while a local maximum appears for 
zero effective field, when B — Bo- In the vicinity of So there 
is a single dimensionless parameter S which determines the 
leading field and temperature dependencies of the thermal 
conductivity. 5 depends on the field via the momentum Ak 
[Eq. 1 37 1] : by substituting Ak into S [Eq. ( 60 1] we obtain an 
approximate (for u 2 ~ vp) expression for 5: 

_ v 2 Ak _ v 2 B eff B e ff 



T 



(70) 



These features can be compared with the effects seen in 
chainJ^ (by interchanging B e j / and B) where the single min- 
imum is at a field B ~ T and the maximum at B — 0. Our 
results for the thermal conductivity isotherms (Fig. |4j display 
similar field and temperature dependance to those measured 
in the experiment!^. 

It should be emphasized that our model relies on some 
simplifying assumptions, and most importantly focuses on a 
purely ID system corresponding to a single ladder. To ac- 
count for the per pen dicular magnetothermal effects measured 
in the experimentf^, our model should be extended to include 
phonons traveling perpendicular to the chains direction. Tak- 
ing into account the coupling of such phonons with the spin 
ladders could result in hybrid spinon-phonon degrees of free- 
dom with higher-dimensional dynamics. Hence, due to this 
hybridization we expect to obtain a higher dimensional spin- 
liquid-like state with strong anisotropics which will account 
for the perpendicular magnetothermal transport. 

An additional limitation on the applicability of our theory 
to a realistic system is that we have assumed a naive model for 
the disorder, and in particular treat it perturbatively. This 
approximation breaks down at sufficiently low T: the disorder 
being a relevant perturbation eventually leads to localization, 
and an effective breaki ng of the ladders to weakly coupled 
segments of finite lengtrpSEO. 

Finally, it should be noted that we have implemented an ap- 
proximate mapping of a ladder onto a chairpS which amounts 
to the truncation of high energy triplet states, and is formally 
justified for J± 3> J\\ ■ Coupling to the high energy sector is 
likely to induce asymmetry between positive and negative de- 
viations of B from Bo , as indeed observed in the experimenlp^'. 

As a concluding remark, in this work we focused on the 
limit vf <fC v, compatible with the parameters of BPCB. 
However, in other quasi ID spin compounds, where J ~ Ob 
(e.g., NaV 2 5 [24] and NO[Cu(N0 3 ) 3 ] [25]), the spinons and 
phonons velocities are comparable in size vf ~ v. Hence a 
strong hybridization between the two degrees of freedom is ex- 
pected in such compounds. Our theoretical approach can be 
extended to account for this phenomenon as well; we expect 
to investigate it further in future work. 
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Appendix A: Umklapp Memory Matrix 

Before proceeding into the calculations of correlation func- 
tions, we show that due to the conservation law (421, simple 



relations between the umklapp matrix elements can be found. 
Substituting Eqs. p7| , p8| into |44f and using Eq. (53]l, we 

get 

Ji + -h = + Ju (Al) 

where Ju is the longitudinal phonons current Ju = 
J dxPidxUi. In addition we have 



Fu = [Ju, H u ] = , =► Fi + F 2 = F^ 



(A2) 



Substituting this into the conservation law ( 42 1 , we find 



relations: 



Fi 


Ak „ 

+ F 2 = — F.. 


(A3) 


see, 


from M pq ~ (F p ; 


F q ), the following 


AO 


Ak 

= - m: 2 , 






Ak 

= =pM s " 2 - M 2 " 2) 


(A4) 


MS 


Ak 

= ^ACi - Ml 2 . 





According to Eq. (51 1 and (52 1, we need to calculate the 



Fourier transform of retarded correlation functions of the form 
C; q (x,t) = </ p u (:M); /, u (0,0)>o, (A5) 
with the force density operators fp(x, t) defined so that 

F; = i[J p ,H u ] = J dxf?{x), (A6) 



in which H u is the umklapp term defined in Eq. (37 1. The 
expectation value (...)o is evaluated with respect to Ho ( 36 1 . 
The first umklapp term to calculate is M" 3 ; from commutator 
identities we find 



f:(x) = i[n(x),H u }= (A7) 
i 9u f d x' m x),co Sm x>) - Akx']] = ^,^,1 - Akx] 

To calculate correlation functions between trigonometric 
functions, we use the result (appendix C in [T]): 



A 3<t>( r i) 



(A8) 



with Ai some constants, K the LL parameter, and T\ = 

\ ln [(^) 2 ( sinh2 (ff ) + sin2 {if))] ' this correlation Ac- 
tion has the property that for J]\ Ai 7^ 0, it equals zero. Since 
Ho is separable in terms of the eigenmodes fa (i — 1, 2), the 
correlation function C" 3 can be written as a product of two 
correlation functions 
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CUx,t) = G\(x,t)C 2 (x,t) = {e 'P«*l(M)-2«*(0,0)l )(e .[2 M2(lll )-2 W2 (0,0)^ 

I 



(A9) 



where the correlation function of each species of the eigen- sponding LL Hamiltonian. This yields: 
modes is calculated independently with respect to the corre- 



Cs 3 (x,t) = sm{-KK 2 j3 2 ) 



To 



2K 2 p 2 +2K 1 c 



(sinh[7rT(a;/ui - t + it)] sinh[7rT(x/ui + t - ie)])~ Kia x 



(sinh[7rT(a;/i;2 - t + ie)] sinh[7rT(a;/«2 + t - ie)])"^ 2 ' 32 , T = 



V2 

2na ' 



(A10) 



with a = 2C and /3 = 2A S defined in Eq. (41 1, the LL pa- where B(x, y) 



rameters K 1 / 2 are defined in Eq. (33 1. The Fourier transform 
Cf 3 (Ak,u) is evaluated using the approximation Kia 2 ~ 0, 
which simplifies the correlation function into an expression 
that can be calculated straightforward by the integraP^I 

I dtsmh-^ 2 (TvTOe^ = ^5(^/4-^,1-^/2) 

(All) 



r(^)r(y) 



is the Beta function. This yields 



C: s (Ak,ij) = 4g 2 u sm(nK 2 p 2 ) 



(A12) 



and consequently the matrix element 



M: a = -i8^ s \ u=0 = J^sm(irK 2 /3 2 )x 

- rp \ 2K 2 l3 2 +2K ia 2 -3 

j Dfr/ /o „T i r^" p[r/ £>2 /n i „•£ i ts Q^^aUT,i-i ^ ^2 / lTf / a 2 , 



T J 
_ v 2 Ak 
= 4ttT 



B[K 2 p 2 /2 -i8,l- K 2 ^]B[K 2 p z /2 +iS,l- #2/9 ]&[*(! - #2/3/2 - i«5) - * (#2/3 /2 + i<5)], (A13) 



After deriving the expression for AC, we wish to show that Using c^e 2 ^* 2 ^'* 5 = 2i/3d x cj>2(x,t)e 2il34 ' 2{x ' t) and integrating 
M" 2 is proportional to M" s , and the rest of the elements are by parts we have 



found from Eq. ( A4 1 . Again, using commutator identities we 
have, 

C: 2 = -iC 1 (x,t)d x C 2 (x,t), (A14) 

then 



M" 2 = / tdt / dxQ[Cs 2 (x,t)]e l 



(A15) 



J dxd x C 2 (x,t)e iAkx = - J dxC 2 {x,t)d x {e iAkx ) = -iAk J dxC 2 e zAkx 



(A16) 
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which gives the simple relation 
1 



M s " 2 

In a similar way 



Ak 



4i (-tAk)Aa => m: 2 = — m: s 



M22 = (^f) 2 M a " s 



(A17) 



(A18) 



By substituting Eqs. (A17l, and (A18l into Eq. ( A4 1 we 



see that all the elements with p or q — 1 vanish, and the 
derivation of the umklapp memory matrix is complete. 



zero. Since we are interested only in the leading temperature 
and field dependencies of the memory matrix, the results of 
the integrals in this section will be important only to get the 
powers of T in each matrix element. The force operators is 
derived from H d [Eq. ([38])] , using 



(Bl) 



Fd = i\J P ,Hd] = i j dx((x)[J p ,cos(2(j>(x))] 



Appendix B: Disorder Memory Matrix 



Now we turn to the calculation of the disorder part of the 
memory matrix. Note that all the non diagonal elements are 



Eq. ( A8 1 is again useful: after disor der averaging, and using 



the identity C(^)C(0) = D8(x) [Eq. (39 1], we get 



J dte^l sinh(TrTt) sinh(— 7rTf ) | ~ 



C S3 (uj) = — 2 



D_ ( T 
^ \2To 



K 2 /3 2 /2+K ia 2 /2 



K 2 /3 2 /2-K ia 2 /2 



D ( T \ K 2/3 2 /2+if 1Q 2 /2-l 



2tt 3 \To 



B[K 2 l3 2 



2nT 



1 - K 2 p 2 /2] 



r 



(B2) 



M d _ .0C SS . n T K 2 p 2 /2+K ia 2 /2-2 

QUI 



D sa 



D 

TTV2 



(B3) 



After some algebra, 



Cti(u)=D 



K 2 f) £ /2+K ia " 12 



dx8(x) / die !Ut x 



\2T J 

^[(sinh[7rT(x/wi - t + ie)] sinh[7rT(»/ui + t - je)])~ Xlct2/4 ](sinh[7rr(x/t; 2 - t + ie)] sinh[7rT(a;/'U 2 + t - ie)])" if2 ' 32/4 , (B4) 

I 



DK-iq 2 ( T x K2fi ' 2+Kia / 2+1 
~~2 {% 



B[K 2 f3 2 /4 + K ia 2 /4 + 1 - — , -1 - K 2 /3 2 /2 - K l0 ?/2] 

4:711 



DK l0 ? 



£>n Si 



2vi 



K 2 fi 2 /2+K ia 2 12 



(B5) 



The result for M 22 is pretty much the same 



dui 
DK 2 P 2 
2v 2 



T \ K 2 p 2 /2+K ia 2 12 



(B6) 
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Appendix C: Phonon-Spinon Memory Matrix ble for phonon-spinon scattering. 

In this appendix we detail the calculation of the correlation 
function appearing in the memory matrix elements responsi- 



C p - 3 (oj) = J dte iwt {F p - s {t);F p - s (Q)) = 

v 4 g 2 p _ s Jdt e ilu+Aw)t S(k-k'+p-p' 

kk' pp' G n ;ki k'^pip^Gni 

K c I' c fc c fc' 1 c *i>( 6 p' 6 p 6 p' 1 6 i>i> - (4i c *i c fc' c *)( & p' 1 & pi 6 p' 6 J>}]- 

Aid = € k — C k - q + Ldp — Ldp + q . 

Using Wick's theorem, one obtains 

{c\,Ckc\, c kl ) = fkfk 1 8 k 'kS k ' lkl ~fkfk 1 S k ' kl S k ' ik +fk 1 S k 'k 1 8 k ' lk , 

(C2) 

and similarly for the other expectation values. We thus get 
C p _ s (w) = ^ W pg 5(w + Au)n p+ qf k - q (l + n p )(l - f k ), 

k P q 

W pq ~ -2v 2 g p _ 3 \p(p + q)\q 2 , q = p - p. (C3) 



G n )5{k 1 - k[ +pi - pi - G nl )pp'(p -p)pip'i(p'i -pi)x 



(CI) 

I 

where / fe = ( e ( e <=-f)/ T + l)" 1 , n p = (e"p /T - l)" 1 , are Fermi 
and Bose distributions respectively. 
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